Magnetized strangelets at finite temperature 
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The main properties of strangelets, namely, their energy per baryon, radius and electric charge, 
are studied in the unpaired magnetized strange quark matter (MSQM) and paired magnetized color 
flavor locked (MCFL) phases. Temperature effects are taken into account in order to study their 
stability compared to the 56 Fe isotope and nonmagnetized strangelets within the framework of the 
MIT bag model. We conclude that the presence of a magnetic field tends to stabilize more the 
strangelets, even when temperature is considered. It is also shown that MCFL strangelets are 
more stable than ordinary MSQM strangelets for typical gap values of the order of 0(100) MeV. 
A distinctive feature in the detection of strangelets either in cosmic rays or in heavy-ion collider 
experiments could be their electric charge. We find that the electric charge is modified in the presence 
of the magnetic field, leading to higher (lower) charge values for MSQM (MCFL) strangelets, when 
compared to the nonmagnetized case. 
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I. INTRODUCTION 

The idea of strange quark matter (SQM) as the ground 
state of nuclear matter, proposed by Bodmer-Witten- 
Terazawa [J-Ojj]; and enforced by Farhi and Jaffe has 
continued alive during the last three decades. Even more 
exciting are the results about the theoretical possibility 
to have the SQM in a neutral paired phase - the so-called 
color flavor locked (CFL) phase [5|,|6j], where all quark fla- 
vors and colors are paired. The true ground state of mat- 
ter would be in this paired phase |6], which is expected 
to occur at high densities and low temperatures. 

If the SQM or CFL phases are the fundamental states 
of matter, one could expect them to be present inside 
compact objects, such as quark or hybrid stars @, [1], 
or in the even more exotic small quark lumps named 
strangelets. The formation of these lumps could be 
caused by the collision of compact stars, supernova ex- 
plosions or primary cosmic rays from strange stars. They 
could even be produced in ultra-relativistic heavy-ion col- 
liders, such as the RHIC 0, and in the Large Hadron 
Collider (LHC) experiments, confirming with their sig- 
natures the formation of the quark-gluon plasma (QGP). 
Furthermore, experiments pointing to cosmic rays are on 
the search for strangelets fiot fllj . 

In early works [J, ll2Hl4| , and in more recent pa- 
pers jTsMlTT l . several studies on SQM in the framework of 
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the phenomenological MIT bag model, and on the prop- 
erties of strange stars and strangelets, have been tackled. 
In particular, they have analyzed the finite size effects of 
strangelets and discussed their stability at T = with 
an energy per baryon lower than the corresponding one 
of the iron isotope 56 Fe, the most stable nucleus known 
in Nature. At T ^ 0, strangelets would be in metastablc 
states, but even if they survive 0.01 milliseconds, their 
signatures would be very important for relativistic heavy- 
ion collider experiments (16J. 

The above studies also suggest that the radii R of 
spherical strangelets have a dependence with the baryon 
number A of the form R = tqA 1 ^ 13 1, where r$ is the 
reduced strangelet radius, which is approximately equal 
to 1.12 fm for the nucleus in the atmosphere. These 
estimates could be of great interest in the search of 
strangelets upcoming from cosmic rays and reaching the 
Earth's atmosphere. For example, the mean free path of 
strangelets in the atmosphere strongly depends on their 
radii [IjJ . 

An interesting characteristic of strangelets arises from 
their relatively low baryon numbers. For A < 10 7 , 
electrons (positrons) cannot be trapped and mixed with 
quarks inside strangelets due to their small Compton's 
wavelength compared to the radius of strangelets [l6| . 
This implies that electric charge neutrality is in general 
no longer ensured, opposite to what happens in strange 
stars. A net nonzero electric charge for strangelets is an 
important parameter to be measured in order to clar- 
ify whether their existence and stability are possible in 
the Universe, and to establish definite conclusions on the 
QGP formation. The sign of the electric charge will also 
give us information about the particle composition, and 
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some authors even attempt to propose possible epic cat- 
aclysmic scenarios depending on it: positively charged 
strangelets would repel ordinary nuclei, but negatively 
charged ones would attract them, turning into strange 
nuclei, so that the normal matter we know would change, 
leading to a great disaster [l6j |. 

Most of the experiments searching for strangelets are 
sensitive only to low baryon number values (A < 10 4 ). 
Thus, it is important to know their theoretical proper- 
ties. Since QCD equations cannot be solved analytically, 
nor by perturbative methods, we shall use in our study 
the bag model, which treats quarks (anti-quarks) as a 
nonintcracting particles confined into a spherical well. 
The latter may be thought as an external parameter 
that realizes the color confinement of quarks and glu- 
ons. Several results about the stability, radius and elec- 
tric charge of strangelets can be obtained from the liq- 
uid drop model (LDM) @, El, widely used in hadron 
physics. In this case, strangelets are described as non- 
interacting Fermi and gluon gases with a modified state 
density, using the multiple reflection expansion method 
(MREM) developed by Balian and Bloch (see also 
Refs. EElMll)- 

In astrophysical scenarios as well as in heavy ion col- 
liders, the magnetic field is an important physical quan- 
tity. Pulsars, magnetars, neutron stars, and the emis- 
sion of intense sources of X-rays could be associated to 
compact objects with an intense magnetic field at the 
surface around 10 13 — 10 15 G [23|, |24|. In the interior 
of these compact objects, the equipartition theorem pre- 
dicts magnetic fields of the order of 10 19 - 10 20 G [25j| . 
Recent studies also indicate that very intense magnetic 
fields, reaching values of 10 19 — 10 20 G, might be gener- 
ated in heavy ion colliders (26[. 

The purpose of the present paper is to study the main 
properties of strangelets at finite temperature in the pres- 
ence of a strong magnetic field. Our analysis is done 
within the framework of the bag model, and using the 
LDM approximation [l5l - fl7l [2l| . By including tempera- 
ture effects, our results can be applicable to astrophysical 
as well as heavy-ion collider environments. In order to 
take account of the finite size effects, we also consider 
nonzero quark masses and the effects of the magnetic 
field and temperature on the surface and curvature en- 
ergy density corrections to the quark and gluon gas. We 
then aim at quantifying the relevant parameters in the 
characterization of strangelets and their stability. 

The paper is organized as follows. In Sec. [TT| we 
present the general thermodynamical expressions to de- 
scribe strangelets at finite density and temperature in 
the presence of a strong magnetic field. We also dis- 
cuss the criterium of mechanical equilibrium based on 
the minimization of the free energy, and briefly address 
the issue of the Debye screening of the strangelet elec- 
tric charge. This effect, although not relevant for paired 
MCFL (CFL) strangelets, it is of major importance in 
determining the total charge of unpaired MSQM (SQM) 
strangelets. In Sec. IIIII we present the numerical results 



for the main properties of MSQM and MCFL strangelets, 
namely, their energy per baryon, radius and charge. We 
also study how the critical baryon number, below which 
strangelets are unstable, depends on several input pa- 
rameters. Finally, our concluding remarks are given in 
Sec. ED 



II. MAGNETIZED STRANGELETS 

AND /i / 

The phenomenological MIT bag model has been wide! 
used to describe strangelets with low baryon numbers ^ 
[T^ - [ll . HU . It has been shown that surface and curvature 
phenomena are indeed important and should be consid- 
ered. In particular, surface and curvature corrections 
to the thermodynamical potential, energy and particle 
density play an important role on their stability at low 
baryon number [16]. At high baryon number, these quan- 
tities approach asymptotically the bulk limit, in which 
the surface and curvature properties are negligible. This 
leads to an unbound phase of SQM, which has been 
widely studied in order to get a more comprehensive un- 
derstanding of the interior of neutron and hybrid stars. 

To consider the effects of the magnetic field on the 
thermodynamical properties of strangelets, we first need 
to know the energy spectrum of the particles involved 
in the system. We assume that our system is under the 
action of a constant and homogeneous magnetic field B, 
guided in the z-direction. The energy levels are quantized 
by Landau levels in the plane perpendicular to the field 
so that the energy spectrum is [27| 



pj 
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where 



P±f 




2n-T] + l) 



B) = 
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to/ are the quark masses, / = (u,d,s), the quantity n 
indexes the Landau level, r\ = ±1 are the quark spin 
projections onto the magnetic field direction, B'j are the 
critical magnetic fields and g/ denote the quark electric 
charges. 

In the absence of a magnetic field, strangelets can be 
considered as spherical drops, characterized by their vol- 
ume V = 4irR 3 /3, surface area S = AirR 2 and extrinsic 
curvature C = 8irR, where R is the radius of the sphere. 
On the other hand, in the presence of strong magnetic 
fields (^ 10 18 G), strangelets are expected to loose their 
sphericity due to the breaking of the spatial rotational 
symmetry, which also leads to the splitting of the parallel 
and perpendicular (to the field direction) pressure com- 
ponents [HI, H3, HI] . Nevertheless, in our study we shall 
assume a spherical shape for strangelets, which turns out 
to be a good approximation for B < 5 x 10 18 G. Deformed 
nonmagnetized strangelets have been considered, e.g., in 
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Ref. [20|| , showing that spherical strangelets exhibit an 
energy per baryon lower than the deformed ones, and 
thus being more stable. 

When dealing with massive u, d, and s quarks, as well 
as with a gluon gas, the general expression for the ther- 
modynamical potential can be written as 



CI 



g t wqq, 

Clg — Clg V V + Clg C C, 

hag v + J2lVf,vV + n Ls s + n Lc c] 
f 



Qqq = Bh 



(3) 
(4) 

(5) 



The quantity fl g takes into account the gluon contribu- 
tion to the thermodynamical potential, while Cl q -q corre- 
sponds to the quark q (anti-quark q) contribution. Glu- 
ons, being massless particles, do not contribute to the 
surface corrections 

HE! SI- up to one-loop approxi- 
mation, the bulk and curvature gluon terms can be ana- 
lytically obtained and read 



dglf_ 

90 



d, 



9 T 2 



36 



(6) 



where d g = 16 is the gluon statistical weight. 

The term J3bag V accounts for the bag energy, with 
-Bbag being interpreted as a vacuum pressure that realizes 
the color confinement of quarks. The sum in the last 
three terms of Eq. ([5]) is over all quark flavors. 

In the one-loop approximation, the bulk statistical con- 
tribution to the quark thermodynamical potential given 
in Eq. §5§ reads 



djT 
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where 



£ = 1 



,-(.E P ,fT»f)/T 



(7) 



(8) 



represents the distribution function of particles (fp) and 
anti-particles (fp) at the temperature T; df — 3 is the 
quark statistical weight and fif is the chemical potential 
per gas flavor. 

The bulk thermodynamical properties of SQM (3^, HH 
and MSQM [13, Hl| have been studied in the astrophys- 
ical context, which is of great importance in describing 
the equation of state of the quark gas, as well as the 
mass-radius relation for quark stars. 

The surface correction effects for quarks given in 
Eq. ([5]) can be written as [HI 



G s — arctan (m//|p|) . 



(9) 



The factor G s takes into account the modification of the 
state density within the MREM. For massless particles, 
it vanishes even at finite temperature. For massive par- 
ticles, their densities on the surface remains always neg- 
ative, i.e., particles shy away from the surface as a conse- 
quence of the boundary condition in the bag model. We 



also note that, in the limit T = 0, and neglecting the 
u and d quark masses, Eq. (EQ) leads to the expression 
obtained by Berger and Jaffe [13j . 

In an analogous fashion, the curvature energy density 
correction to the thermodynamical potential of quarks 
can be written as [la, uSSl : 



f,c 



G r =1 



G e ln(f+f j 
arctan (nif /\p\) 



(10) 



48tt 3 

In this case, the factor G c has been derived [l6[ to fit the 
shell model. 

Since the anisotropy in the spectrum, imposed by the 
presence of a magnetic field, quantizes the energy states 
of every particle, the integration over dp x dp y in Eqs. ([7]). 
© and (fTO)) are to be replaced by the rule 



oo r+oo 



oo J — oo 



dp x dp v -> 2nq f B > C 11 ) 

J7=±l n=0 



where the sum over Landau levels is limited up to 



I 



x f ~ 1 
2B/B c f 



(12) 



due to the fact that the Fermi momenta must be real- 
valued quantities [13] ■ Here we have denoted by Xf — 
Hf/mf the dimensionless quark chemical potentials and 
I[x] the integer part of the real number x. 

The role of the magnetic field in the quantities associ- 
ated to surface and curvature phenomena has been stud- 
ied in Ref. [34j. Due to the presence of a magnetic field 
these quantities exhibit divergences for the ground-state 
Landau level (n = 0). Indeed, when p±j = 0, the sur- 
face contribution of Eq. (|9]) diverges logarithmically as 
^2fClf. s ~ — nln(p± u p±dP±s), while the curvature en- 
ergy density given in Eq. (TP0"j) diverges as Ylf^f,c ~ 
-fi(l/pxu + l/pxd + 1/pxs) + ^\n(p ±u p ±d p ±s ). Since 
these divergences only appear in the limit p±f — 0, they 
can be removed by introducing an infrared cutoff in p±/LJ. 
Above this cutoff, the surface and curvature terms are 
then numerically well-behaved quantities. 

To study the stability of strangelets the free energy 
must be minimized. This minimum represents the most 
stable configuration that the quark droplet can reach. 
The free energy can be written as 



F 



F„; 



Fg — Qg.v V + flg iC C, 

Fqq = B bag V + £ [F f>v V + F f>s S + F fjC C] 

f 



(13) 



1 We note that the inclusion of the anomalous magnetic moment 
of quarks in the energy spectrum acts as a natural infrared cutoff 
f°rpx/[12- 
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where 



F f:V = Qf, v + fifNf, v , 
F f,s = Qf, s + fJ>fN fiS , 
F f,c = fi/, c + HfN f , c . 



(14) 



Bulk, surface and curvature corrections to the particle 
density are 



N 



N 



on 



dfif ' dfif J ' u dfj,f 

respectively, so that the total particle density reads 

Nf = N Lv V + N Ls S + N Lc C. (16) 
In turn, the baryon number A is 

^=~E^/- 



N 



(15) 



(17) 
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Equilibrium configurations can be found by solving the 
equation 



dF_ 
W 



= 0. 



(18) 



N,T,B 



The energy of the gas is obtained from the free energy 
by adding the entropy contribution, i.e., E = F + T S. 
The entropy S contains the bulk, surface and curvature 
contributions, 



S f,v 



f,v 



dT 



on 



dT 



s 



f,c 



dT 



(19) 



respectively. 

As electrons (positrons) are forbidden to coexist with 
quarks in strangelets due to their small radii, electric 
charge neutrality is never ensured. This is an important 
feature in their detection. The electric charge, derived 
from the free charge distribution, is given by the expres- 
sion Z[ ICC = z v + z s , where 



E Qf N f,» ' 

f 



£ 



1fN L cC. 



(20) 
(21) 



For strangelets, however, the effects of the Debye 
screening are not negligible and should be taken into ac- 
count [14, 35, 36]. The strangelet charge distribution will 
depend on the Debye screening length A^> , 



\ D 2 = 4nJ2l 2 f 



dN 



(22) 



By solving Poisson's equation in the Thomas-Fermi ap- 
proximation, the screened volumetric charge is then cal- 
culated as (l4| 
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Z v = — RX 



D u 



1 — — — tanh f 
R \Xd 



(23) 



and the total electric charge becomes 

Z = Z v + z s . (24) 

The corresponding Coulomb energy of the strangelet is 
given by 



E c = ^ 2 z 2 v X A D 



3X D , / R\ , _ 2 / R 
—— tanh - — + cosh - — 
R \Xd J \Xd 



(25) 



which should be added to the total energy of the system. 

Finally, it is worthwhile to mention that, in the study 
of strange stars, one can derive from Eq. (J7J all the ther- 
modynamical bulk properties and the equation of state 
(EoS) of the system. In this case, the baryon number A 
is a huge quantity (A ~ 10 57 ) so that the surface and 
curvature contributions can be safely neglected. On the 
other hand, to describe the mechanical equilibrium con- 
figurations of strangelets, where A < 10 7 , the finite size 
effects, like the surface tension in Eq. ^ and the curva- 
ture energy density in Eq. (fTTj|) . play an important role 
in determining their properties. These terms are also 
important in the study of the quark-hadron phase tran- 
sition gi m us . 



III. NUMERICAL RESULTS: MSQM 
STRANGELETS AND MCFL STRANGELETS 



In this section we solve numerically Eq. (|18p . which 
guarantees the mechanical equilibrium of strangelets. 
Our aim is to study the energy per baryon for strangelets 
in the MSQM and MCFL phases as a function of the 
baryon number A, bag energy density the magnetic field 
B, temperature T and gap parameter. We shall com- 
pare our results with the case of zero magnetic field, and 
also compute the radius of strangelets and their electric 
charge as a function of the baryon number. First we con- 
sider strangelets formed from unpaired MSQM, and then 
we study those corresponding to the paired MCFL phase. 



A. Unpaired phase: MSQM strangelets 

Strangelets with A < 10 7 require zero electron and 
positron number densities. If they are formed from SQM 
or MSQM, the /3-equilibrium condition implies the equal- 
ity of chemical potentials: fi u = = fi s . In the presence 
of a magnetic field, however, the latter condition does 
not imply the equality of bulk particle densities Nf lV in 
MSQM strangelets. Indeed, we recall that, at T = 0, 
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FIG. 1. (Color online) Energy per baryon for unpaired MSQM 
and SQM strangelets with B 



75 MeV fm J and consider T - 



5 x 10 1B G. We take B bs , e 
-■ 0, 15 MeV. 



FIG. 2. (Color online) Critical baryon number A c[it of MSQM 
strangelets for 23 = 5 x 10 18 G and T = 0, 15 MeV, as a 
function of Bbag- The nonmagnetized SQM case is also shown. 



correspond to the z-component of the magnetic Fermi 
momenta and the magnetic mass, respectively [27j- The 
presence of hJj in the Fermi momenta and B c ^ in the parti- 
cle densities implies that, despite the chemical potential 
equalities, the particle densities should be different for 
every flavor. 

Equation (|18p can be written in the form 



dF, 
dV 



dF n -, 



N,T,B 



8V 



0. 



N.T.B 
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dV 
dF al! 



dV 



N.T.B 



N,T,B 
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(28) 



(29) 
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dE, 



and 



bag 



c 



dV 



R~ 



R 2 



fi c , (30) 



(31) 



The above system of equations, together with Eq. dTTl) . 
allows us to determine the common quark chemical po- 
tential and the radius of the strangelet, as well as to 
evaluate other quantities such as the energy per baryon, 
particle densities, and electric charge. 

In Fig. [TJ we present the energy per baryon for un- 
paired MSQM strangelets at zero and finite tempera- 
ture, assuming m u = rrid = 5 MeV, m s — 150 MeV, 
B bag = 75 MeV fin" 3 and B = 5 x 10 18 G. The en- 
ergy per baryon of SQM (£> = 0) is also shown for com- 
parison. We observe that MSQM strangelets exhibit an 
E/A lower than SQM strangelets. The horizontal dotted 
line corresponds to the iron energy per baryon. We as- 
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FIG. 3. (Color online) Radii of strangelets as a function of the 
baryon number for B b a g = 75 MeV fm" 3 and T = 0, 15 MeV. 



Radii of MSQM (B = 
strangelets are depicted. 



5 x 10 18 G) and SQM (B = 0) 



is a consistent a ppr oximation in the magnetic field range 
B< 10 19 G [II, III. 

The presence of the magnetic field reduces E/A even 
at T = 15 MeV, for the bag parameter and quark masses 
considered. Finite lumps of magnetized quark matter 
with high baryon numbers could be absolutely stable 
at this temperature, which opens the possibility of find- 
ing strangelets from the coalescence of quark or hybrid 
stars where those high magnetic fields are expected to be 
present. Clearly, temperature effects increase E/A due 
to the increment of the thermal energy of quarks and 
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FIG. 4. (Color online) Ratio Z/A 1/3 (in units of e) for B = 
5 x 10 18 G and T = 0, 15 MeV. We take S bag = 75 MeV 
fm -3 . The nonmagnetized SQM strangelet electric charge is 
also shown. 



gluons. 

For both unpaired MSQM and SQM strangelets, the 
energy per baryon as a function of the baryon number 
can be fitted with the general expression [161 ] 



E E s E c 



(32) 



where the coefficients E v , E s and E c are ascribed to 
bulk, surface and curvature phenomena, respectively. 
This fitting function shows that, for high baryon num- 
bers, the energy per baryon approaches its bulk value, 
E/A ~ 3[if+TS v /A. For low baryon numbers, the terms 
depending on A show destabilizing properties due to their 
growth [HI- As expected, the bulk limit increases 

with temperature. Surface and curvature terms get bal- 
anced with each other to maintain the free energy at its 
minimum for a fixed A. 

We also observe from Fig. Q] that, in the presence of a 
magnetic field, strangelets can be stable for a wider range 
of the baryon number. Furthermore, depending on the 
bag parameter, temperature, and magnetic field, there 
exists a critical value A CI i t with E/A CI n = 930 MeV, such 
that for A > A cl -i t strangelets can be in a stable phase and 
long lasting in Nature, but metastable when A < A cv - lt . 
In Fig. [21 the critical baryon number is plotted as a func- 
tion of the bag parameter for B = and B = 5 x 10 18 G, 
and T — 0,15 MeV. It can be seen that for sufficiently 
large -Bbag, there are no stable strangelets. For B = 0, 
stability requires B hag < 70 MeV fm~ 3 at T = and 
Sbag £ 66 MeV fm- 3 at T = 15 MeV. For B = 5x 10 18 G, 
the upper limits are Bb ag ^ 77 MeV fm -3 at T = 0, and 
Sbag < 73 MeV fur 3 at T = 15 MeV. The allowed range 
of bag values is thus wider in the presence of the mag- 
netic field. It becomes clear that the magnetic field plays 



a crucial role on the stability of strangelets. This is ob- 
served in bulk, surface, and curvature terms. On the 
other hand, temperature effects turn out to be less sig- 
nificant in comparison with those of the magnetic field 
on the surface and curvature terms. Notice also that, if 
E/A > m n ~ 939 MeV, SQM strangelets would exhibit 
neutron emission [l3l . [To ] . This opens the interesting pos- 
sibility to detect strangelets by the emission of particles 
from unknown astrophysical sources. 

The dependence of the radii of SQM and MSQM 
strangelets on the baryon number is depicted in Fig. [3] 
Due to the thermal energy of quarks and gluons, temper- 
ature effects tend to slightly increase the radii. The same 
behavior is observed in the presence of a magnetic field. 
This is a direct consequence of the decreasing behavior of 
the surface and curvature terms as discussed above. The 
energy, preserving the stability on the surface, is dimin- 
ished by the magnetic field, which tends to stabilize more 
strangelets, but also increases their size compared to the 
nonmagnetized ones. As a function of the baryon num- 
ber, the strangelet radius R behaves as R = tqA 1 / 3 , with 
ro ~ 0.9 fm, for high values of A. Temperature and mag- 
netic field tend to modify the values of tq in about 1%, 
being for MSQM greater than the corresponding ones of 
SQM, but always smaller than the nuclear value. 

Finally, the dependence of the total electric charge of 
strangelets on the baryon number is depicted in Fig. [4] 
The electric charge increases with A and, for a fixed A, it 
also increases with the magnetic field. For small baryon 
numbers, Z oc A, while for large A the relation between 
Z and A differs from the linear one. This is caused by the 
screening of the bulk's electric charge for strangelets with 
a radius greater than the Debye length. The strangelet 
electric charge is distributed within a husk of a Debye 
length radius from the surface. The relation Z oc A 1 / 3 is 
observed [3], when the contributions of the free surface 
electric charge are not considered. Including these effects, 
as in Eq. (|2"4")) . we obtain the asymptotic behavior Z oc 
A 2 ' 3 . 



B. Paired phase: MCFL strangelets 

The stability, EoS and mass-radius configurations of 
a quark gas in the MCFL phase have been studied in 
Ref. [40( within the framework of the bag model. In anal- 
ogy to unpaired SQM and MSQM strangelets, we now 
analyze the relevant properties of strangelets composed 
by matter in the MCFL paired phase at finite tempera- 
ture. In the limit B = and T = 0, our results agree 
with those obtained by Madsen in Ref. 



22] 



We assume that quarks behave like noninteracting par- 
ticles, but the binding energy of the diquark condensate 
is associated to a BCS pairing energy gap A. Following 
Ref. , we take a common value of the gap parameter 
for the predominant color pairings (ud,us, ds) [4ll |42||. 
Furthermore, in order to simplify our analysis, we ne- 
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FIG. 5. (Color online) Energy per baryon for CFL [B = 0) 
and MCFL (B = 5 x 10 18 G) strangelets at T = 0, 15 MeV. 
We fix Stag = 75 MeV fm" 3 and A = 100 MeV. 



FIG. 7. (Color online) Critical baryon number A C iit of CFL 
strangelets as a function of the strange quark mass m s . The 
plots are given for 23 = and 5 x 10 18 G, T = and 15 MeV, 
taking Stag = 75 MeV fm" 3 and A = 100 MeV. 




FIG. 6. (Color online) Critical baryon number A C rit of CFL 
strangelets as a function of Bbag- The plots are given for 
B = and 5 x 10 18 G, T = and 15 MeV, and assuming 
A = 100 MeV. 



FIG. 8. (Color online) Critical baryon number Actit of CFL 
strangelets as a function of the gap parameter A. The plots 
are given for B = and 5 x 10 18 G, T = and 15 MeV, 
assuming Bbag = 75 MeV fm -3 . 



gleet any dependence of A on the magnetic fielc0. On 
the other hand, as our study is devoted to strangelets at 
finite temperature, we consider a temperature-dependent 
gap parameter in analogy to superconductivity studies of 



quark matter [2l|, |44|, |4£ 



A = 2~ 1/3 A 



1/2 



(33) 



2 A more accurate study would require the determination of A by 
solving the corresponding gap equations which depend on the 
magnetic field [3Sll43|. 



where T c = 2 1 / 3 e 7 Ao/7r ~ 0.71A is the critical tem- 
perature of the superconducting phase, above which the 
pairing of quarks is forbidden. 

For the magnetized CFL phase, the bulk thermody- 
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FIG. 9. (Color online) Radii of nonmagnetized and mag- 
netized CFL strangelets for a magnetic field value B = 
5 x 10 18 G, T = 0, 15 MeV, B b a g = 75 MeV far 3 and 
A = 100 MeV. 



FIG. 10. (Color online) Ratio Z/A 2/s (in units of e) of 
nonmagnetized and magnetized CFL strangelets for B = 
5 x 10 18 G, T = 0, 15 MeV, B bag = 75 MeV fin -3 and 
A = 100 MeV. 



namical potential in Eq. (|31l) reads as 

3AVb 



^ - 1! 



f 



(34) 



where the second term in the right hand side accounts 
for the cost of the free energy which is compensated by 
the pairing formation [i(J, and \xb = (/■*« + Md + Ms)/3 is 
the baryon chemical potential. 

This color symmetric state requires quark bulk particle 
numbers to be equal to each other to minimize the free 
energy and enforce the color and electric charge neutral- 
ity of the CFL phase [f| I45T - I48} . This requirement leads 
to the equations 



fJ-B 



N, 



2A 2 MS 



N, 



'Mb 



(35) 



From Eqs. ([33)) it is clear that the bulk electric charge 
for CFL and MCFL strangelets is zero, and the only con- 
tribution comes from the surface. Charge screening is 
thus negligible for CFL strangelets. This is in contrast 
with the charge behavior of SQM and MSQM strangelets, 
for which the bulk charge contribution is dominant. 

In an analogous fashion to the unpaired phase, the 
free energy density, which now contains the quark pairing 
contribution, has to be minimized by solving the equilib- 
rium equation ([TBI) , or equivalently, Eqs. (|28)) - (|3"Tj) . with 
Ec = and Q, v given by Eq. (IM)) . In Fig.[S]we present the 
energy per baryon E/A as a function of the baryon num- 
ber, for nonmagnetized and magnetized CFL strangelets. 
The results are given for B = 5 x 10 18 G, T = 0, 15 MeV, 
B bag = 75 MeV far 3 and A = 100 MeV. We notice 
that the energy per baryon, corresponding to strangelets 
in the paired phase, has a lower value than for unpaired 



strangelets. Furthermore, MCFL strangelets exhibit an 
E/A lower than CFL strangelets. 

As in the case of unpaired quark matter, surface 
and curvature terms are responsible for the stability of 
strangelets at low baryon numbers, and the energy per 
baryon can be fitted with a general expression of the 
form given in Eq. ([32]). At high values of A, surface 
and curvature phenomena become negligible and E/A 
tends to its bulk value. Furthermore, bulk contribution 
increases with temperature, while surface and curvature 
get balanced with each other to maintain the free en- 
ergy at a minimum. For given values of the bag and 
gap parameters, temperature, strange quark mass and 
magnetic field, the ratio E/A of CFL strangelets always 
fixes an A cr \t above which the stability is guaranteed 
(compared to E/A — 930 MeV). Figure [5] shows the de- 
pendence of A cr ; t on the bag parameter for B = and 
5 x 10 18 G, T = and 15 MeV, and taking A = 100 MeV. 
Once again, for sufficiently large -Bbag> there are no sta- 
ble strangelets. For B — 0, stability requires B^ ag < 
105 MeV far 3 at T = and B hag < 102 MeV far 3 at 
T = 15 MeV. For B = 5 x 10 18 G, the upper bounds are 
Sbag ~ 112 MeV fm ~ 3 at T = 0, and B hag < 110 MeV 
fnr 3 at T = 15 MeV. Therefore the presence of the mag- 
netic field leads to a larger stability window. 

To see the sensitivity of A cr n to the variation of other 
parameters, we present in Figs. [7] and [8] how the critical 
baryon number varies with the strange quark mass m s 
and the gap parameter A. From Fig. [7| it is clear that 
A cr it is a slowly increasing function of m s . On the other 
hand, Fig. [8] shows that j4 cr i t is a decreasing function of 
the gap A. In both cases, the effect of the temperature 
is to increase the required values of -A cr it in comparison 
with those at zero temperature. 
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The CFL and MCFL strangelet radii at T = 0, 15 MeV 
are shown as a function of the baryon number in Fig. [51 
for fixed bag and gap parameters. It is observed the same 
behavior than that of the radii of the unpaired strangelets 
studied in the previous subsection, i.e. R ~ 0.9A 1 / 3 . 
Furthermore, the effect of the magnetic field is to yield a 
larger strangelet configuration due to the contribution of 
surface and curvature effects. This is in contrast with 
magnetized compact stars, for which surface and cur- 
vature corrections are negligible, and the magnetic field 
yields star configurations with smaller radii. We also note 
that an increase of A would lead to higher values of the 
radii, as in the case of compact objects [40(. 

Finally, in Fig. [10] we present the dependence of the 
electric charge on the baryon number for CFL and MCFL 
strangelets. The curves show that temperature always 
tends to decrease the electric charge. In the bulk, the 
equality of quark densities implies zero electric charge 
contribution. Therefore, these strangelets are character- 
ized only by their surface electric charge distribution. As 
can be seen from the figure, the magnetic field also tends 
to decrease the surface electric charge. When A is large, 
the behavior of Z as a function of A is Z ~ 0.3A 2 / 3 [22[ 
for nonmagnetized CFL strangelets, while its character- 
istic value is Z ~ 0.2A 2 / 3 for MCFL strangelets. 

IV. CONCLUSIONS 

Strangelets from MSQM and MCFL have been stud- 
ied within the MIT bag model and the LDM framework. 
A comparison with SQM and CFL strangelets for mas- 
sive quarks has been presented. We have examined how 
the stability and size of quark lumps are modified by 
the temperature and the presence of a strong magnetic 
field. While temperature tends to destabilize strangelets, 
the magnetic field acts in the opposite direction. Any 
nonvanishing temperature contributes to an increase of 
the energy per baryon E/A, which measures how sta- 
ble are strangelets when compared to the isotope 56 Fe, 
or any other given pattern, such as neutrons, pions or A- 
particles. On the other hand, the magnetic field enhances 
their stability, both in bulk and surface. Strangelets in 
the MCFL and CFL phases show more stability than un- 
paired strangelets, as expected from the gap energy and 
pairing of quarks. 

The stability of strangelets at zero and finite temper- 
ature strongly depends on several parameters and, par- 
ticularly, on the bag constant, strange quark mass, gap 
parameter, and the magnetic field. There exists a crit- 
ical baryon number A CV i t below which strangelets are 
metastable states, and can decay and radiate particles 



from their surface, even at T = 0. The critical baryon 
number increases with B^ ag , m s and T, but is a decreas- 
ing function of the magnetic field B and the gap parame- 
ter A. These results, nevertheless, should be taken with 
care due to the approximate nature of the liquid drop 
and bag models. 

The radii of strangelets were also studied in order to 
have an estimate of their size. As the pairing gap in- 
creases, the radii of strangelets get larger, as in the case 
of compact objects [4(| • Magnetized strangelets turn out 
to have a radius bigger than nonmagnetized ones, since 
the magnetic field tends to relax the surface tension, thus 
increasing the strangelet radius. 

Finally, we have also found that the charge neutral- 
ity of strangelets is never ensured and, in the presence 
of the magnetic field, strangelets are always positively 
charged. Therefore a nonvanishing electric charge may 
constitute the main characteristic to be searched for in 
experiments or in the detection of strangelets from astro- 
physical sources. Unlike MCFL (CFL) strangelets, which 
are neutral in bulk, unpaired MSQM (SQM) strangelets 
have a contribution to the electric charge due to the 
screened bulk charge. The relation (Z, A) represents an- 
other characteristic that could distinguish the phase of 
strange quark matter present in strangelets. For large 
A, MSQM and SQM strangelet electric charge exhibits 
the behavior Z oc A 1 / 3 (attaining Z oc A 2 ! 3 asymptoti- 
cally), whereas for MCFL and CFL strangelets, the ratio 
Z/A 2 / 3 reaches a constant value. These could be striking 
properties in their detection either in cosmic rays or in 
heavy-ion collider experiments. 
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